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o 

CsJ ' Eigenvalues of the Breit Equation 

U ' 

< 

QQ ' in which only the static Coulomb potential is considered, have been found. Here the 

^SJ , detailed discussion on the simple caces, ^Sq, m = M and m y^ M is given deriving the 

exact energy eigenvalues. The a^ expansion is used to find radial wave functions. The 
' (-H ' , leading term is given by classical Coulomb wave function. The technique used here can 

O ' be applied to other cases. 



(dip + l3im)aa'Spi3' + 5aa'{-a2P + /32M)/3/3' Saa'SjSfS' 



*„,^, = E^ap, 



i^H , Subject Index The eigenvalues of the Breit equation 

^: 

1. Introduction and the Breit equation 

^ , 1.1. Introduction 

The Breit equation [1] has traditionally [2] considered to be singular and has never been 



in 

'sj" , tried to solve itself. Instead, the Pauli approximation or a generalized Foldy-Wouthysen 

^ ' transformation has been applied to derive the effective Hamiltonian H^, consisting of the 



Schrodinger Hamiltonian 



^° = -2U^mJ^ "7 

t,j . and many other relativistic correction terms [2] [3]. Then perturbation method vv'as used to 

^ , evaluate the eigen-values (binding energies) in power series in a as high as possible (or as 

- - -' much as needed for experimental verification of QED). 

Here we shall try to directly solve the Breit equation. The radial wave fundtion for the 
simplest case, ^S leptonium states with equal masses, is given by the confluent Heun function. 
The auther is not familiar to handle this function. Instead, the approximate method using 
a^ ~ 10~^ expansion shall be used to find wave funcitons. Its leading term is given by 
Coulomb wave function. 

However, the exact eigenvalues are found in singlet-states with equal m = M and unequal 
masses m ^ M cases. We shall study the simplest case, ^S, equal lepton masses, m = M, in 
greater details here. The technique used here shall be applied to other cases, ^{l)i, ^{l)i and 
^Pq including m ^ M, to find exact eigenvalues. 
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More complicated cases, m ^ M described in §2.3 below and triplet cases, shall be 
described in due course. 



1.2. Basic equation 

The Breit equation for two leptons (Dirac particles), with charge and mass, (— e, m) and 
(+e, M), interacting through the (static) Coulomb potential — e^/r = — a/r (a=l/137) is 
given by 



{a.ip + I3im)fj,u5pa + Sfj,u{-a2P + ^M)^ 



-SiLuSn 



^,. 



E^,p. 



(1) 



This equation holds in the CM system. The total energy of the system E can be expressed 
as 



E = ^J' 



m^ 



+ v/m^ 



The momentum operator in (1) can be written as 



^ 1^ 1 9 

P=-V = -— , 
I I or 



(2) 



(3) 



introducing the operator r canonically conjugate to p. Then the (static) Coulomb potential 
(contribution form one longitudinal photon exchange between two oppositely charged point 
Dirac particles) can be written as — e^/r = —a/r. aj,l3j{j = 1,2) are the usual 4x4 Dirac 
matrices for particle 1 (— e,m) and 2 (+e,M), respectively. ^ai3 is 4 x 4 Dirac spinor wave 
function. 
We shall try to solve 



(ai-V + /3im) + (-52- V + /SaM) - - 
t I r 



^ = £;^, 



and find the eigen-value(s) 



E= ^Jm^ -q^ + \fM^ 



(4) 



(2) 



which is specified by discrete real value(s) qn- 

The two particle system can be classified by the total angular momentum and its z- 
component and parity. In so doing it is sufficient to specify spin angular parts of the large- 
large components ^ a^ (a, /? = 1>2) to be 

\l)h (Z = 0,l,2,--.) 

3(/-l)z +3 (/ + !),, (/ = 1,2,...) 
and 3(/),, (/ = 1,2,---) 

where I is the quantum number for the orbital angular momentum. Spin-angular parts of 
other (large-small, small-large, small-small) components of '^ap are completely fixed by those 
of the large-large components. 
We shall describe the singlet case in §2 in detail. Other cases can be handled similarly. 
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2. The Singlet Cases 

2.1. ^{l)i, general 

We shall first discuss the simplest case of ^{l)i. ^afs can be taken as follows: 



a 


/3 


*a/3 


1 
2 


1 
2 


F(r)|i(0^> 


1 
2 


3 

4 


.{G(r)|3(/ + 1),>+G(r)|3(/ -!),>} 


3 

4 


1 
2 


i{/7(r)|3(Z + l),>+#(r)|3(Z -!),>} 


3 

4 


3 

4 


K(r)|i(0« > 



(5) 



where \^{l)i > and p(/ it 1)/ > are normalized spin-angular wave functions (z-component of 
the total angular momentum has been omitted). The Breit equation (4) in §1 gives the fol- 
lowing set of differential equations for the radial wave functions F{r),K{r), G{r),G{r),H{r) 
and H{r): 



a 

E^ m- M] F 

r 



^a-^K..-) 



I 



21 + 1 \dr 



d l-l 



(G-H), 



^i(^y-±^)iG-H) + VTTT(^^-^-^')iG-H).o, 



E + --(7n + M)] F 
r 



a 

E-\ ^m- M] H 

r 

a 



E+- + (m + M) 
r 



l + l d I 



21 + 1 \dr r 



K, 



(F + K) 



E^ m + M]G, 



a 



E-\ \-7n- M] H 

y r / 

a 



I 



21 + 1 \dr 



"^'^lUF + K) 



E + --m + M]G. 
r 



(6) 



Therefore, we have the relations: 



E + - + m-M)H=(E + --m + M]G, 
r J \ r 

E + - + m-M)H=(E + --m + M]G. 
r J \ r 
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For later convenience, we introduce the dimensionless quantities: 

P = 2gr, 



E ^M2 - g2 + y^ 



g- 



2aq 



2aq 
_ M + m + E 

2aq 
_ M + m- E 
2aq 



(>0), 
(>0), 
(>0). 



(7) 



If the eigen- value i? is a proper relativistic generalization of the Schrodinger case (as we 
shall see it is the case), y and A are large quantities of the order of 1/a^, while u is the order 
of unity. 

For p — )• cxD(r — )• cxd) all radial wave function should behave like 

F(p)^p"{l + 0(l/p)}e-^/2. (8) 

This is because e^f'"^ = e"^"^ and E = ^Jm? — q^ + -^/M^ — q^. Notice that if we consider 
E, the total CM energy ^Jm? + p^ + y M^ + p^^ to be larger than M + m, we would expect 
that the radial wave functions should contain a factor e^^ . If we analitically continue p to 
imaginary values we get the factor e"'^'' and the energy expression (2) for bound states. 
The singlet case becomes very simple when two masses are equal M = m. First 

H = Cr, H = Gr, 



K 



1 — up 



F 



Then, the combination 



obeys the differential equation 



h" + h'l-l + - + 



l + Xp 
1 1 



(9) 



(10) 



p ' p{'^ + yp)j 

\\2 J p^ 4 p^ 


l{l + l) 

p2 


1 


2p{l + yp)\ 



+h 



2.2. ^Sq, M = m 

The ^5*0, i.e., I = 0, case we have the equation for h 

p p(i + yp) 



+h 



a^y 



1 



p 4 p2 '2p{l + yp)_ 



0. 



fill 



If we ignore a^ and higher order terms in (11) and notice that a^y is of the order of unity, 



we find 



h" + h' 



P 



a^y 



11^=0, 
P 



in) 
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which is precisely the Schrodinger equation for our issue and h in this approximation is given 
by F{—^^ + 1,2, p) = F{1 — n,2,p) apart from the normahzation constant. ^^ — 1 must 
be equal to an integer n — 1 = 0, 1, 2, • • • . 

(11) shows that p = is the regular singular point while p = oo is an irregular singular 
point. 

Near p = h{p) can be expressed as 

h{p) = p'h{p), 

, 1 (13) 



where s = —1 + a/ 1 . 



Another solution s = —1 — \ 1 — ^ is unacceptable from the square integrability of the 
wave function ^ai3- h{p) obeys the equation: 



p p{l + yp) 

h{p) can be expanded in a Taylor series near p = 0: 

oo 

h{p) = Y,Kp'', (15) 

n=0 

I { a^y 1 1 , 

(n + l)(?i + 3 + 2s)-/i„+i 

y 

( 1 a^v 1 1 

+ {n{n + l + 2s)--{n -^ + 1 + s + - + sy)} K 

[ y 2 2 J 

n-l-^ + l + s\ K-i = 0. (16) 

It can be shown that h{p) can not be a polynomial: 

TV 
h{p)=Y,hnp'' {N< 



N 

ooj 

n=0 



from the recurrence formulae (16). So the sum in equation (15) must extend to n = oo, a 
sharp difference from the Schrodinger case. Note that the expansion (15) holds only the tiny 
region p < 1/y ~ 0{a'^). 
Next we discuss on p = oo, which is an irregular singular point. Assuming the form 



h{p) = e^V' 



Co H \ h 

P P2 



at p — ;• OO (17) 



where A, k, Cn are constant. Introducing (17) into the differential equation (14), we find A is 
either or 1 while k is undermined. From the square-integrability of the wave function we 
must choose A to be and k < oo. 
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k turns out to be 



k 



a^y 



1-s 



(18) 



appeared in (14), whose value is not fixed from the situation at p — t- oo. Further discussion 
will be given in §2.2.4. 

2.2.1. Detailed discussion on ^Sq, m = M. Note again 
h{p) obeys the diffrential equation, (14), 



7 1 . \ N 
-l + -\+h 



d^h dh 
dp'^ dp 



where 7 = 2 + 2s, s = -1 + Jl - ^ 



1 



p(i + y/3) \dp 



d 



h, 



4 ' 



S = h + sy, 



N = ^-l-s, 



p = 2qr, 



(19) 



^ 2Vm^-g^ ^ _E_ 
" 2aq 2aq ' 

a^y and sy are of the order of unity 0(1). 

The RHS of equation (19) is of the order of 1/y or O(a^) ~ 10"'^. As the first step we shall 
ignore the RHS of (19). Then h can immediately solved: 

/i^F(-iV,7,p). (20) 

Normalizability of the radial wave function requires N equals to positive integer: 

cP'v 
N = — ^-l-s = n-l, n = l,2,3,---. (21) 

n is the pricipal quantum number. (20) is the relativistic extension form of the non- 
relativisitic Coulomb wave function 



N = ^-l-s ^n-1, 

F(-iV,7,p) ■ 

7 = 2 + 2s ^2. 

The energy eigenvalues derived from the condition (21) are 



(22) 



Er, 



2\J w? — q^ = 2mi / 1 



1 



1 + {ayy 



2m Jl 



a^ 



{o^yY + OL^ 



2m Jl 



a^ 



{2(n + s)}2 + a2' 



(23) 
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When n = 1, 



o? 



E, = 2m^l--. (24) 

The approximate solution F{—N,j,p) is correct ignoring terms of the order 0{l/y) ~ 
a^ ~ 10~^. We write 

hip) = Fi-N,^,p) + fip). (25) 

Whenn = 1, F{0,-/,p) = 1, 

/(O) = 0. (26) 

We shall introduce 

h = F + f, 

into the eq. (19), then we find 



^ + ^1-1 + 21 + /^ = / ^fi 

dp'^ dp { p] P \ pC^ + yp) \dp 



^ + ^A-^ + -}+f- = \-7T^rz::,{i:-^]}(F + fl (27) 



whose solution gives h = F + f. 

Or, alternatively, we may use an approximate method: F is the order of unity, while f{p) 
is the order of 0{l/y) ~ O(a^), so that / in the RHS of (27) can be ignored. This equation 
gives immediately the solution 



f{p) 



7o (Pyp^ J [ a{l + ya)\ da 



6}F 



It is difficult to perform integretion here, but computer shall easily do the jobs. 

The equation (21) give the exact eigenvalues for our problem. The reason is as follows. 
When p —7- CO, the function h{p) = F{p) + f{p) behaves like (see §2.2.4 ) 

hip) = p^(co + ^ + ---Y P = N 

as will be described in §2.2.4 F{p) contributed to p" , while f{p) does not, i.e., 

— 3- ~ O ( — I , when p — )• oo. 

pf^ \ypj 

Therefore (21) gives the exact eigenvalues. 
To find /(p), we may use 1/y expansion, or a^ expansion, since y is of the order of 1/a^. 
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Introducing 

oo 

/(0)(p) = F{p) = F{-N,r,p), N = n-1, 

oo 






where f^ip) is of the order of {l/yY ~ 0{a ), and the differential operators 

dp^ V Pj dp p 



where D^'^' is of the order of {^/yY (here z^ = 0, and 1), the differential equation (19) will 
give the following serries of equations 

Z?W/W = 0, 

D(0)/('^)(p) = Z)«/(-"i)(p), (z. = l,2,...). 

It is easy to see that 

/(°)(p) + /«(p) = F(p) + /«(p) 

shall be good if the terms of the order of a^ ~ 10~^ are ignored. 
These techniques can be applied to other singlet and triplet cases with m ^ M. 

2.2.2. f{p) for the ground state of^So, m = M, near p = (p < l/y). The function h{p) 
for n = 1 is given by 



hip) =F{0,j,p) + fip) = l + fip), 

/(O) =0. 
The Taylor expansion of f{p) near p = is given by 

oo 

/(p) = E/"/'" 

n=l 

/n's are equal to /i„,(n > 1) and given by (16). Noting that 



sy = -l, S = h 



(28) 



one finds 



N = 0, 7 = 2, 



/o = 0, /i = Y^7 f2=- ^yfi = -^y> 



8/14 



and in general 



In 



+1 



-yr/2 



(n + l)(n + 2)(n + 3)' 



fip) ^ E/"/'" = E 



i-y) 



n~\ 



n=\ 



^^n(n + l)(n + 2)2' 

n=l 



2y [ 2(yp)2 



log 



1 \ 3 ^ 
1 + ypy 4 2yy9 



which holds at — - < p < -. f{p) is of the order of 0{l/y) = O(a^) and this yf{p) is good 
within the errors of 0{l/y) = 0{a'^). General cases n > 1 can be treated analogously. 

2.2.3. h{p) for the ground state of^So, m = M, near p = —1/y, mathematical curiosity. 
The eq. (14) has the regular singular point at p = — 1/y. The regular solution near p = —1/y 
is given by 



h{x) = ^ gn.i 



n=2 

where x = p — -. 

Another solution near x = = p contains log of the solution (29). 

The recurrence formulae for g^s are 

3g3 = {2(l + 7)y + 2-7}52, 



(29) 



(30) 



-{n + l)(n - l)gn+i + {n{n + 'yy)6} gn-y{n-l- N)gn-i. 
Ignoring 0{\/y) terms, one finds 



(31) 



Taking 92 = 1, 



53 = 2y52, 9a = ^yg3 = 3y^52- 



9n+l t 7 7T5n " ^y" \ {u > 1). 



(n-1)- 



So that h{p) turn out to be 



hip) = ^gnx^y^-'=^nx-+'y^ 

n=2 n=l 

x'^y [l + ypf 



1 — yx y'^p 

This expression is valid at the unphysical range —2/y < p < and good within a^. 
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2.2.4. Behavier of h{p) at p ^ 00. p = 00 is the irreguler singular point, so that h{p) can 
be expressed in the form (17) as described in §2.2. 
We shall give the recurrence formulae here for the case oi^S,m = M 

h{p) = /fco + ^ + ^ + -- 

ci = |-/3(/3 + 7-l) + J}co, 

(n+l)yc„,+i + (n-/5){(n + 1 -7-/3)y-n + (5}c„ 
+ (n-l-/3)(n-7-l-/3)c„_i = 0, 

where 

7 = 2 + 2s, 

S = i, 

2' 

2 

B = N = -^-l-s = n-l. 
^ 2 

N has already been fixed in Eq. (21). 

We repeat that fJ^ is the contribution from F{p) and f{p) = h{p) — F{p) does not 
contribute to p^ term. Furthermore it is easy to show 

U \ — when p — )• 00. 



p^ \yp 

This shows (21) is in fact the correct (or exact) eigenvalues. 

2.2.5. Comparison of Breit case with Dirac Coulomb solution. For simplicity we shall 
discuss only ground states of ^5*1/2 (Dirac) and ^Sq (Breit). 
Dirac equation with Coulomb potential [4] reads 



|iaV + /3m- fW = E"^, 



(32) 



E= Vrn2 
^ for ^-51/21 spin up, can be written 



^ ^ / F{r)Yoo^^ \ si/2 spin up 
I K{r){af)^ J pi/2 spin up. 



(33) 



0-3 



1 





, Yoo{e,<P) = 


1 





V27r 




ar = 


-(ar). 
r 





The differential equation for F is given by 

d^F ( 1 a 2] dF f/^ a\2 2 1 r. ^ 

+ \-Emr—---o + -)-r- + \(E + -] -m^ F = 0. (34) 



dr'^ \E + ^ +m r^ r j dr 
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Introducing dimensionless quantities 



P = '^qr, y 



E + m 
2aq ' 



and write F to be 



F{p) = e"P/'p^h{p), 



where s is introduced to remove 1/p^ terms as usual. Then the differential equation for h{p) 
is given by 



," + v;-i + i±^+ ' 



p p{l + yp) 



1 



+ '^^T-'-V-p-[^^''jpii + yp) 



0, 



(35) 



where E = \Jvii^ — g^, binding energy B = m — \Jvrfi — q"^, y = {E + m)/2aq. However s is 
different from Breit case 



-1 + yr 



Q^ 



The form (35) in the Dirac case looks exactly the same as the Breit case (14). Only 
differences are in below. 





Dirac 


Breit 


s 

N 






aE 

1 — s 



However the ground state solution of (35) in the Dirac case is 

h{p) =F{0,2 + 2s,p) = l, for n = 1. 
And both 



aE \ ,1 

1 — •? =0 and — \- sy = 

q J 2 



holds and they give the same value of q, 



q = am. 



Therefore 



E = V^2(l_Q,2) 



B = m 



E = [l-Vl-a^} 



m. 



(36) 



(37) 



As is well-known, Dirac equation with Coulomb force can be solved for all cases of n and 
I in terms of two confluent hypergeomotric functions (multiplied by p^e'P'"^). 
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However one cannot do 



A^ = and 5 = 



simultaneously in the Breit case. 



a^y 



then 



While 6 = demands 



A'^ = means 1 — s = 0. 

2 






1 a^ 

7^ y = 

2 8 ^ ■ 



but 



2 ■ 



2 1st exiced state! 



and N cannot be zero. 

From these situations, the Breit case is more complicated than the Dirac case, as described 
in this article. 

2.2.6. Confluent Heun function. The solution h{p) of differential equation (14) is conflu- 
ent Heun function[5], which has 5 parameters. Unfortunately the author is not familier to 
this function, so that the auther has used physical ground and 1/y expansion of the wave 
function h{p). 

2.3. ^Sq, M ^m 

Unequal mass case M ^ m can be treated in exactly the same way as in §2.2 for the case of 
M = m. 
The equation for the quantity h{p): 



F + K = h{p)e 



-p/2 



is given by 

h" + h' 



-1 + 



+ 



2{1 + yp) 



p p{l+yp) p{(l +yp)2 - (Af -m)2p2} 



+h 
+h 



+ 



l + yp 



p 2p{l+yp) p{(l+y/>)2-(M-m)2p2}J 



1 a\^ 
4 + T^^ 



(M + m)2p2 ) [ (1 + yp) 

(i + yp)2 



P' 



(M - friY 



0, (38) 



where M = M/2aq and fh = m/2aq. 

As in §2.2, p = and p = oo in eq. (38) are the regular and irregular singular points, 
respectively. Therefore, we can adopt the same procedure as in §2.2. 
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To remove the l//0^ terms, put 



h{p) = p'h{p), 



finding again s = — 1 + ^Jl — cP^/A. The differential equation for h{p) reads 

1 { {m- Mfp 1 + yp 



(A) 
(D) 



h" +h' 
+h 
+ 



-1 + ^ + 



p (1 + ypY -(fh- Mfp^ \ l + yp 



+ 



1 



(1 + ypY - {fh- Myp^ \ l + yp 



+ 



P 



1 



{fh-Mfp l + yp 



2 [ (1 + ypf -{m- Mfp^ \ l + yp 
+\ + '^{y'-2{rh^ + M')] 



+ 



P 



+ 



o? im^ - M2)2 



-P 



0. 



4 {l + ypf 

We should notice here that (A)+(B) and (C) reduce 

-8 



1 



(^) + (^) ^ W^^r '=2+^^' 



p 



{m — MYp l + yp\\ s 



P P'' 



(39) 



(C) -^ 0, and [D) = 0, 

when the equal masses case, M = m. (39) reduces to (14) for the case, ^5o,rra = M. 
We shall rewrite the term (D) as follows: 



(D) 



a^ (m^ - M^y ( yp 



[E)' 



(F)... 



4 y^ 

4 I y2 



l + yp 



yp 



l + yp J p 



r2^x2 



a'^y ( {m^ - M'')X I 
'~2 



y" J p 

The last term [F) shall be added to the first line of th coefficient h[p) in (19): 



h" + h'{-l + ^-2{B)\ 



+ h 

+ h[ 


'o?y\^ (fr?-M''\^\ ^ 


- s 


1 


(A) + [B) + (C) + {E)\ = 0. 


P 



(40) 

(A), (B), (C) and (E) are of the order of Ijy or o? . Now we can use the same technique 
described in §2.3, ^S, m = M. 



hip) = F(-iV,7,p) + /(p), 



(41) 
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N 
N 



a^y \ ( fr? -M^^ ^ 



2 [ ^ ^^ 

n — 1, n = 1,2, • 



1-s, 



(42) 



n is the principal quantun number. The background part / is 0(o?^ and can be found by 
solving 

/" + /'{-l + ^}+/(^) = |+2(i?)^-(A)-(i?)-(C)-(£)|F(-iV,7,p). 

Thus the accuracy of the solution F + / is quite good, terms of the order of a^ are ignored. 
Other radial wave functions, i^, K^ G, -ff , G and K can now be writen down. 
The energy eigenvalues from (42) 



o?y 



1 



fr? - M 



2\2^ 



aE 
1^ 



1 



n? -M 



2\2" 



n^ s ^n^ 



E = ^W^ - g2 + ^m^ - q2 
is easily calculated. Therefor the binding energy B = M + m — E \s 



B = M + m- E= ^ 



m=l 



a \ 



2m 



(43) 



ei 



62 



63 



2mM 
(M + m) 

2mM 
(M + m) 

4mM 



3mM 
{M + m)2 



1 



5mM 



+ 



f){mMf 



(44) 



(M + m) 1^^ (M + m)2 (M + m)"^ 

and so on. The first term is exactly the well-known form of the reduced mass. 
It is easy to see that (43) reduces to (23) when M = m. 



Knowing F + K = e ^''^p^h{p), all relevant radial wave functions, F,K,G,H,G and H 
can easily be derived. 
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